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Abstract 

We analyze the higher-twist effects and the S'C/(3)-flavour symmetry breaking in 
the correlation functions used to calculate form factors of pseudoscalar mesons in 
the QCD light-cone sum rule approach. It is shown that the Ward identities for 
these correlation functions yield relations between twist-4 two- and three-particle 
distribution amplitudes. In addition to the relations already obtained from the QCD 
equations of motions, we have found a new one. With the help of these relations, 
the twist-4 contribution to the light-cone sum rule for the pion electromagnetic form 
factor is reduced to a very simple form. Simultaneously, we correct a sign error in 
the earlier calculation. The updated light-cone sum rule prediction for the pion form 
factor at intermediate momentum transfers is compared with the recent Jefferson 
Lab data. Furthermore, from the correlation functions with strange-quark currents 
the kaon electromagnetic form factor and the K — > ir weak transition form factors 
are predicted with 0{m s ) ~ 0{m 2 K ) accuracy. 
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1 Introduction 



An accurate knowledge of the pion and kaon light-cone distribution amplitudes (DA) in- 
troduced in the studies of hard exclusive processes in QCD ][]] is important for various 
frameworks where these DA are being used. Among the most topical applications one 
could mention the calculations of exclusive semileptonic and hadronic 5-meson transitions 
into pions and kaons using pQCD 0, QCD factorization [[J or light-cone sum rules [|J. 
Although there are definite indications that at the normalization scale of O(lGeV) the 
leading twist 2 pion DA is already quite close to its asymptotic shape, one still encoun- 
ters a large uncertainty of the nonasymptotic part. Moreover, very little is known about 
nonasymptotic SU (3)-flavour asymmetry in the twist 2 kaon DA. 

One of the promising ways to study DA is to employ vacuum-to-pion or vacuum-to- 
kaon correlation functions of light-quark currents. At high virtualities, using the operator- 
product expansion (OPE) near the light-cone, these correlation functions are expressed in 
terms of DA. On the other hand, the same correlation functions are related, via dispersion 
relations, to the observable form factors of pions and kaons with the contributions of 
excited hadronic states approximated by quark-hadron duality. In the resulting relations, 
known as light-cone sum rules (LCSR) ||, the experimental data on form factors can be 
used to yield nontrivial constraints on DA. The LCSR for the pion electromagnetic (e.m.) 
form factor was derived in Refs. || [7| and for the 7*771"° form factor in Ref. [Q. In order 
to further increase the accuracy of these sum rules one has to gain a better control over 
higher-twist effects in the OPE. In the case of the pion form factors the most important 
subleading contribution to the LCSR is of twist 4. The kaon e.m. form factor which so 
far was not analyzed in the LCSR framework, demands also inclusion of the twist 3 effects 
proportional to m s . 

The aim of this paper is twofold. First, we analyse the higher-twist effects in the 
vacuum-to-pion and vacuum-to-kaon correlation functions. We demonstrate that a new 
useful tool is provided by standard Ward identities for the conserved e.m. and axial (in the 
chiral limit) currents. Simultaneously, we correct a sign error in the previous calculation of 
the twist 4 term and update the LCSR prediction for the pion e.m. form factor. Second, 
we include S77(3)-flavour symmetry breaking effects at 0(m s ) ~ 0(m 2 K ) in the correlation 
functions. We calculate the twist 3 part and obtain LCSR for the kaon e.m. and K — > 7r 
weak transition form factors at intermediate spacelike momentum transfers. 

The plan of the paper is as follows. In Sect. |2| we introduce a generic correlation 
function, which yields LCSR for the pion, kaon, and K —>■ ir form factors for different 
flavour combinations of light-quark currents. The correlation function is then calculated 
with twist 4 accuracy including first-order in quark mass terms. In Sect. |3] we derive the 
Ward identities in the chiral limit and demonstrate that they lead to relations between 
two- and three-particle DA of twist 4. In Sect. |] the numerical results for the pion form 
factor are presented with a corrected twist 4 contribution. A comparison of our prediction 
is made with the recent data on the pion e.m. form factor obtained at CEBAF. Sect. [5] 
contains LCSR results for the kaon electromagnetic form factor, and Sect. ^ deals with 
the K — > 7T weak transition form factor. We summarize our conclusions in Sect. [F| The 
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Figure 1: Diagrams corresponding (a) to the leading- order of the correlation function 
(b) to the contributions of twist 3,4 quark- antiquark-gluon DA. Solid, dashed, wavy 
lines and ovals represent quarks, gluons, external currents and pseudoscalar meson DA, 
respectively. 

appendices contain the expansion of the quark propagator in App. [A], the definitions of the 
DA and their asymptotic expansions in App. |B] and in App. |Cj the as corrections to the 
twist 2 LCSR obtained in Ref. 0. 

2 Correlation functions 

As a starting point, we introduce a generic correlation function: 

T„v(p,q)=i J rf 4 a;e^(0|T{(g 2 (0)7 M 75gi(0)) {e l q l {x) lv q' 1 {x)+ e 2 ^(x) 7 „g 2 (a;))}|P(p)) , (1) 

where, in order to obtain the LCSR for the pion e.m. form factor the following quark- 
flavour combination has to be taken: qi — q[ — u, q 2 = q' 2 = d. In this case the on-shell 
hadronic state P = 7r + , and t\ = e u = 2/3, e 2 = e<i = ~ 1/3 are the quark e.m. charges in 
the units of e. To calculate the kaon e.m. form factor, one simply has to replace d — > s and 
7r + — > K + in the above. There are two other, physically interesting correlation functions 
yielding two independent LCSR for the K — > n weak transition form factors obtained from 
Eq. (]]]) at qi = s, q[ = u, q 2 = d, P = n + , e± = 1, e 2 = and at q% = d, q 2 = u, 
q 2 = s, ei = 0, e 2 = 1, P = K°. Summarizing, if one calculates the correlation function 
(0) the result can easily be adjusted to any of the flavour combinations listed above. If the 
external 4- momenta squared q 2 and (p — q) 2 are spacelike and large, the operator product 
in the correlation function (HD can be expanded near the light-cone in terms of pion or 
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kaon DA of increasing twists. One may then retain a few first terms in this expansion, 
having in mind that higher twists are suppressed by inverse powers of Q 2 = —q 2 and/or 
\(p ~ l) 2 \ (f° r a more detailed discussion see e.g. Refs. |7], |||). There are two leading-order 
diagrams obtained from the two terms in Eq. ([[]) by contracting the quark fields q\ with q\ 
and q2 with q 2 , respectively and replacing them by the free-quark propagators. The first 
diagram proportional to t\ is depicted in Fig. [L]a. The second diagram, proportional to e 2 
is obtained from the first one by changing the direction of the quark line and replacing the 
quark-flavour indices 1 <-> 2. The next-to-leading approximation for the quark propagator 
generates the diagram in Fig. |l|b (and its ~ e 2 counterpart) which brings three-particle 
quark- ant iquark-gluon DA of twist 3 and 4 into the game. This diagram is calculated using 
the first-order in gluon field term in the light-cone expansion of the quark propagator given 
in App. HI. We systematically retain all terms of 0(m q ) ~ 0(m 2 p ) in order to be able to 
account for SU(3) breaking effects in the LCSR for the kaon form factors. At the same 
time, the 0(m 2 ) contributions arising, e.g. from the denominators of quark propagators 
are neglected. 

The result for the correlation function (|T|) obtained to twist 4 accuracy reads 
i 

T^{p, q) = ifp J du\Ti(Q 2 , s, u)p^p u + T 2 (Q 2 , s, u)p^q u + T 3 (Q 2 , s, u)q^p u 
o 

+T 4 (Q 2 , s, u)q^q u + T 5 (Q 2 , s, u)g^ (2) 

with 

T x {Q\s,u)^ 2u[eilpp{u) - e ^ p{U)] 



uQ 2 — us + uum 2 P 
1 [4/ 3 p 



J Vai [eim qi (p 3P (ai) - e 2 m q2 (f 3P (ai 



{uQ 2 — us + uurrip) 2 \ ff 
-2u ^4 [ei^ip(u) - e 2 gip(u)} - 4 [e 1 G 2 p(u) - e 2 G 2P (u)} - 2u [eiflf 2 p(«) + e 2 g 2P (u)} 

u 



(1 - 2v) \ 2[eiip ±P (ai) + e 2 <f±p(ai)] + \ei(p\\ P (ai) + e 2 y?||p(a; 
-2[ei<^_Lp(ai) - e 2 (p± P (ai)] - [ei^p^) - e 2 (p\\p(ai)]\) ) , (3) 

T 2 (Q 2 ,s,u) = 



._> , eiV9p(w) - e 2 (f P (u) 



uQ 2 — us + uump 
+ (uQ 2 -us\uum 2 P ) 2 {f ~ 

+4 [eig 1P (u) - e 2 gip(u)] - 4 [eiG 2 p(u) - e 2 G 2P (u)} - Au [e 1 g 2 p{u) + e 2 g 2P (u)] 



+ / Va { 



4(1 - v) [e 1 (p ±P (a i ) + e 2 ip_Lp(ai)] + (1 - 2v) ei<p\\ P (ai) + e 2 ip\\ P (a 
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-4(1 - v) [ei(f± P (ai) - e 2 (p± P (ai)] - e x <p\\ P {ai) - e 2 tp\\ P (ai) 



(4) 



T 3 (Q 2 ,s,u) 
4 



uQ 2 



eiip P { u) - e 2 fp{u) 
us + uump 

\ —ftp 



e^m^tpapiu) - e 2 m q2 (p aP (u)) 



uQ 2 — us + uump) 2 \ 3 
+4 [e x g 1P (u) - e 2 g XP {u)) - 4 [e x G 2P (u) - e 2 G 2P (u)} - Au [e x g 2P (u) + e 2 g 2P (u)] 



+ / Va, 



- Av [e-i^pxpioii) + e 2 <p ±P (ai)} + (1 - 2v) e x ip\\ P {ai) + e 2 (f\\ P (ai) 



-Av [eiV9_Lp(ai) - e 2 ip±p(ai)] 



ei(p\\ P {oci) - e 2 ip\\ P {ai, 



(5) 



T 4 (Q\s,u) =4^ 



e x g 2P (u) + e 2 ^ 2 p(«)] 



2 \2 ' 



T 5 (Q 2 ,s,m) 



■uQ 2 — us + uump) 
Q 2 + s + (u - u)mp 

2T [CiVpW - e 2 ip P (u)\ 



(6) 



2(uQ 2 
fi P 



(uQ 2 — us + uum 2 P ) 
Q 2 + s + (u - u)m 2 p 



us + uum P ) 

- [e x m qi tp pP {u) - e 2 m q2 ip pP (u)} 



2(uQ 2 



us + uump) 2 



4 [e x g XP (u) - e 2 g 1P (u)} - 4 [e x G 2P (u] - e 2 G 2P (u)} 



u 




+ J Voii 


(l-2v) 



eiV||p(ai) + e2^||p(Q!i 



ei^pta*) - e 2 (p\\ P {a i/ 



(7) 



where s = (p — g) ; u = 1 — it; aj = ai, a 2 , 1 — «i — a 2 ; «i = «2, a X} l — a x — a 2 and 

1i M 1— U 



M — «! 



1 — a x — a 2 



1 — «! — a 2 



In the above, y>p is a generic notation for the twist 2 DA of a pseudoscalar meson P = tt 
or whereas y9 p p, <^ CT p, v?3P and g XP , g 2P , <p± p \\ p , <p± p \\ P are, respectively, DA of twist 
3 and 4. Their definitions taken from Ref. [O (see also Ref. O]) are collected in App. [BL 
The decay constant /p of P is defined as (0 | q-zlnlsli I -P(p)) — ifpPn- Furthermore, 
/ip = m 2 p /{m qi + m g2 ) is the twist 3 DA normalization factor and G 2P (u) = Jq dvg 2P {y ) . 
In the case of nonstrange quarks, q x = q' x = u, q 2 = q' 2 = d, both chiral and isospin 
symmetry limits can safely be adopted. In this limit the u, d quark masses as well as the 
pion mass are neglected and DA are either symmetric or antisymmetric (see App. [B]) [], 

lr The type of symmetry is established applying G-parity transformation to the underlying matrix ele- 
ments. 
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with respect to the replacements u u, or a\ «-> «2- In this case the twist 3 parts in 
Eqs. vanish and the combination of quark charges e x — e 2 = e u — e<i = 1 factorizes 

out. The resulting expression for T\ coincides with the one obtained in Ref. || except that 
the signs of the terms containing the twist 4 quark-gluon DA (p±p j \\p are opposite. The 
same discrepancy in signs is found in the expressions for Tj obtained in the chiral limit 
in Ref. |12[ comparing them with Eqs. ©-(0). In the next section we will demonstrate 
that Eqs. ©-(0) are fully consistent with the relations obtained from QCD equations of 
motion. Finally, we note that the twist 3 terms in Eqs. (0)-(0) are new. 

3 Ward Identities 

Multiplying the correlation function ([!]) by the four-momentum q one obtains 

q u T, u = - /rf 4 xe^^(0|T{g 2 (0)7 M 7 5 gi(0)^-(e 1 g 1 (a;)7^(^ + erf 2 (x)jM*))}\P(p)) 

-*(zoK0|[g 2 (0)7 M759 i(0^^ , (8) 

where the second term containing equal-time current commutators originates from the 
differentiation of the 9{xq) in the T-product of currents. For the conserved vector currents 
qi = q[ and g 2 — 9 2 the ^ TS ^ term on the r.h.s. of Eq. (JJ) vanishes. For the second term 
the standard commutation relations for the equal-time current densities can be employed, 
e.g., in the case of the pion: 

^(0)7^75^(0), (e u u(0, x)7 w(0, x) + e d d(0, x)~/ o d(0, x)) 

yielding for the correlation function the Ward identity 

q u T^ = %Uvn- (10) 

In the chiral limit m qi = m q2 = the axial-vector current is also conserved. Hence, we get 
an additional relation 



= 5 (3) (x)d(0,f) 7At 75«(0,f), (9) 



(p-qYT^ = -if nPu . (11) 

The above Ward identities are valid for arbitrary q and p. This circumstance allows one to 
get relations between various pion DA by substituting Eq. @ in l.h.s. of Eqs. (|l0l) ,(|il"l). 

Here we will only concentrate on the chiral limit, so that both Eqs. ([TI]) and ( |TTD 
are valid and p 2 = 0. It is easy to check that the r.h.s. of these equations are saturated 
by the twist-2 contribution to their l.h.s. Hence, the Ward identities ( P~0| ) and fllTD yield 
nontrivial relations between two- and three-particle DA of twist 4. Note that in the chiral 
limit different twists are separated by dimensions, therefore contributions to the correlation 
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function with twist higher than 4 neglected in our calculation are unimportant^. Using 



2q.p d 1 



(q — up) 4 du (q — up) 2 



and 



1 d 
+ u- 



(q — up) 4 (q — up) 2 du (q — up) 2 



together with partial integration in u, rewriting all twist 4 contributions in Eqs. (|i"0|) and 
([□]) as Jq 1 du ^ uq y F(u) and then extracting F(u) = 0, one obtains the following relations: 



92* W 



T>ai((<pxTr((Xi) - (1 - 2v)<f± 7r (a i )) 



(13) 



G 2 Ju) 



u 



Va{o ((p±_ n (ati) + (p± n (ai)) 



(14) 



g l7T {u) = G 2n {u) + -uug' 2n {u) 



-- / V<\; 



(l-2v) (ip\\ n (ai) +2ip ±7T (ai)\ - if^iai) -2ip ±7T (ai) , (15) 



where g' 27C {u) = dg 2 {u)/du. We notice that the above expressions can be used to rewrite 
in the chiral limit (P = n) the twist 4 part of the correlation function (§) using only one 
DA g 27T and its derivative over u, so that 



T^ u = if 7, du\y—^-—{2up^p u - q^p u - p^q v - (q ■ p)g^)yj 



u 



v {2u 2 g 2lT (u) - 2u 2 ug 27r (u)^j + (p^q v + q^p v ){~2ug 2lT (u) + uug 2w (u)) 



(uQ 2 — us) 2 

+ {p^u ~ %Vv + 2g M g v )0hr(tO + fiV [( 2 Q 2 + (Q-P)0- + 2u )) 92k{u) - (q.p) uug' 2% {u)\ \ .(16) 



The relations (|13j) and (|T4| ) can also be obtained using the technique of QCD equations 
of motion pTOf] . The starting objects in this case are the derivatives of quark-antiquark 
operators expressed via quark-antiquark-gluon operators, e.g.: 

^- (0\d{0) lfll5 u{x)\7r + (p)) = i£ada(o\d(0)^ 5 x x G Xt "(ax)u(x)\7r + (p)) (17) 

and 

^- (0\d(0) lfllf3lul5 u(x)\7r + (p)) = % £ada (o\d(0) lf , lf3lul5 x x G Xu (ax)u(x)\iT + (p)) , (18) 



2 In fact, we also neglect four-particle Fock components of twist 4 in the light-cone expansion of the 
matrix elements. This is consistent with the approximation adopted in deriving the relations from QCD 
equations of motion [pi. 



where G^ u = g s G®(\ a /2), tr\ a \b = 25 . The relations derived from Eq. ( [T7| ) are 



and its u <-> u equivalent: 



(19) 



#1tt(«) 



G 2 n(u) 



v) (ip\\ v {oti) - 2(p ±7T (ai 



(20) 



Combining Eqs. (|T9|) and fl20"D one gets the two relations obtained in Ref. [ IQ| . Eq. ( p"8D 
which was also used in Ref. yields 



^lTr(tt) 



m 1 

2#27r(w) + GW(u) + - 



VaiV (y\\ w (ai) + 2^(0^)) 



(21) 



and 



9iM 



u 1 U 

^2tt(w) + G 27r (u) + 2 / 2^(1 ~ v ) (-^lk(«i) + 2<^j 



.a,- 



(22) 



Combining Eq. (|2~T1) and Eq. (p2|) with Eqs. ( |19|) and (|20"1), after simple algebra one indeed 



reproduces the relations obtained from Ward identities, but only two of them, Eqs. ([13]) 
and (|l4l) . The relation (|l^), the only one involving the DA ip\\ n is new and was not 
obtained using equations of motion. Note that the observed consistency between the 
relations derived from Ward identities and from QCD equations of motion provides an 
independent check of our result for the correlation function. Indeed, taking the correlation 
function calculated in Ref. [12] with different signs at the terms with <p± nt \\ n we obtain a 
contradiction between the two types of relations. 

If the chiral symmetry is violated, m q ~ m 2 P ^ 0, the Ward identity (|H]) based on the 
conservation of the e.m. current is still valid, but the relations following from this identity 
are more complicated, mixing DA of twist 2,3 and 4 and not allowing to reduce the twist 
4 part to an integral over a single DA. The corresponding analysis goes beyond the scope 
of this paper. 



4 Updated prediction for the pion e.m. form factor 

The LCSR for the pion e.m. form factor was originally derived in Ref. ||. Let us briefly 
outline the procedure. The part of the correlation function @j (in the chiral limit) propor- 
tional to ~ ViiPv was matched to the hadronic dispersion relation in the variable s = (p—q) 2 , 
that is, in the channel of the axial- vector current: 

■f},T(n* \ 2if rr F 7r {Q 2 ) , J p h (s')ds' 

if n / duTi{Q ,s,u) = h / — ; • (23) 

J —s J s — s 

s h 
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In this relation, the first term on the r.h.s. is the ground-state contribution of the pion 
where f n = 132 MeV and F n (Q 2 ) is the pion e.m. form factor defined in the standard way: 

F n (Q 2 )(2p - q) v = (7r + (p - q) | f u m | rc + (p)) , (24) 

jem \) e [ n g q Uar k e .m. current. The contributions of the a\ meson and excited states 
with J p = 0~, 1 + form the spectral density p h which is estimated as usual, with the help 
of quark-hadron duality: 

. , i 

p h {s)Q{s - sj) = ^ ( dulm s T 1 (Q 2 ,s,u)e(s-s^), (25) 

71 J 


where the effective threshold parameter Sq = 0.7 GeV 2 is determined from the SVZ sum 
rule |13| for the correlator of two u^^^d currents. Representing the l.h.s of Eq. ( ^3|) in a 
form of the dispersion integral: 

1 °° A I 1 

I ' duT 1 {Q 2 ,s,u) = - [—?— fduInvTUQV,^), (26) 

J 71 J S' — S J 



using Eq. ( E2H ) and performing the Borel transformation , (p — q) 2 — > M 2 , we obtain the 
resulting sum rule: 

s o i 

F^{Q 2 ) = —J dse~ s ' M2 J du Im.TUQ 2 , s,u) . (27) 



o 



In the twist-2 approximation one has 0: 

Fl 2) (Q 2 ) = £ju ^(n,/.)exp [~^\ • (28) 

In the above Uq = Q 2 /(sq + Q 2 ) and we have indicated the dependence of the DA ip n on 
the normalization scale /i. 

The LCSR ( p7|) was further improved in Ref. where the 0(a s ) contribution to 
the twist 2 part was calculated by taking into account the perturbative gluon exchanges 
between the quark lines in the diagram of Fig. |l|a. For convenience we present in App. 
the explicit expression for F^ 2 ' as \Q 2 ) which has to be added to the r.h.s. of Eq. (|2"8|). 
Recall that this contribution provides the ~ a s /Q 2 asymptotic behavior |j] of the form 
factor. As explained in detail in Ref. J7| the form factor obtained from LCSR includes both 
the hard-scattering and soft (end-point) contributions. 

Our main update of the sum rule for F n (Q 2 ) concerns the twist 4 term for which a new, 
corrected expression is obtained from Eq. (|T6|): 
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where 

¥^ 4) 0>^) = 2u(g 27T (u,y) - ug' 2n (u, fi)) . (30) 

The second term on the r.h.s. of Eq. fl29|) is a 'surface term' originating after the 
continuum subtraction as explained in Ref. |7|. In the same paper the factorizable twist 6 
contribution to LCSR was calculated: 

Fi6)(g2) = T^ (0| ^ |0)2 (31) 

in terms of the quark condensate density (see Ref. [[jj for the diagrams and other details 
concerning this contribution). Note that the twist 6 term is numerically very small starting 
from Q 2 = 1 GeV 2 which is therefore a natural lower boundary of the LCSR validity 
region Q 

We turn now to the numerical calculation of the pion form factor, 

F n (Q 2 ) = Fi 2 \Q 2 ) + F^\Q 2 ) + F^{Q 2 ) + F^{Q 2 ) , (32) 

where twist 2,4 and the factorizable twist-6 contributions to LCSR are added together. In 
our numerical evaluation of Eq. (|32"1), following Ref. @ we take 0.8 < M 2 < 1.5 GeV 2 and 
adopt the variable normalization scale /z 2 = (1 — u)Q 2 + uM 2 of the light-cone DA. The 
same scale is adopted for the normalization of a s . For the latter the two- loop running is 
used with A( 3 ) = 340 MeV corresponding to a s (l GeV) = 0.48. For the twist 2 pion DA we 
take the asymptotic form <p n (u) = 6m (1 — it). The influence of nonasymptotic corrections 
will be discussed later. Concerning the twist 4 pion DA we actually need only one of them, 
g2n- Interestingly, in first order of the conformal expansion Jl0| this DA does not contain 



nonasymptotic contributions. Using the asymptotic form of ^(w) presented in App. [TJ 
one obtains a compact expression: 

20 

^(u, y) = -5 2 Muu(l - u(7 - 8u)) . (33) 

The nonperturbative parameter <5 2 ~ 0.2 GeV 2 determining the vacuum-to-pion matrix 
element of the quark-antiquark-gluon current (see definition in App. 0) was estimated 
from various 2-point QCD sum rules in Ref. ]0|. To assess the theoretical uncertainty, we 
have recalculated 5 2 using the diagonal sum rule for two quark-antiquark-gluon currents 
which is less dependent on the variations of quark and gluon condensates. The result, in 
agreement with Ref. |f0)| , is 

5 2 (lGeV) = 0.17 ± 0.05 GeV 2 (34) 



obtained with (0 | qq | 0) = (-240 ± 10 MeV) 3 and (0 | {aJ^G^G 11 ^ | 0) = 0.012 ± 
0.006 GeV 4 . 

Our prediction for the pion e.m. form factor given by Eq. ( p2|) is plotted in Fig. ^|, 
calculated with the asymptotic pion DA at a typical value of M 2 = 1 GeV 2 , and at y — y u 
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Figure 2: Pion e.m. form factor obtained from LCSR with the asymptotic pion DA (solid) 
including the twist 2 leading-order (long- dashed), twist 2 0(a s ) (short- dashed), twist 4 
(dash-dotted) and factorizable twist 6 (dotted) contributions. 

and 5 2 = 0.17 GeV 2 . Importantly, the corrected twist 4 contribution is about two times 
larger than estimated before || [7| . Note that at Q 2 — > oo the twist 4 term given by Eq. ( p9|) 
has the same ~ l/Q 4 asymptotic behavior as the twist 2 contribution (^) [J but has one 
extra power of 1/M 2 . This can be seen explicitly by rewriting Eq. (|2"5|), with the help of 
Eq. (p3|), in the form of a dispersion integral with the integration variable s = Q 2 u/u: 

vm-^l^j&vU-S^-i). (35) 







yielding at Q 2 — ► oo 



f<W)~^(l-e-*« 2 ), (36) 



3 Recent work on the pion and kaon form factors at low momentum transfers can be found in Ref. [ fL4| . 

4 Contributions nonvanishing or growing with Q 2 are absent in LCSR. Such anomalous contributions 
emerge in QCD sum rules based on the local condensate expansion (l^, |l7j making the latter not applicable 
at Q 2 > 1 GeV 2 . 
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Figure 3: The pion e.m. form factor calculated from LCSR in comparison with the 
Jefferson Lab data ]7j|/ shown with points ( the experimental error and the model uncertainty 
are added in quadratures) . The solid line corresponds to the asymptotic pion DA, dashed 
lines indicate the estimated overall theoretical uncertainty ; the dash-dotted line is calculated 
with the CZ model jl^ / of the pion DA. 

to be compared with the corresponding limit of Eq. (p8|): 

Fi 2) m ~ ^ (l - (l + J,) e-oV" 2 ) (3?) 

Although the twist 4 term has indeed an extra suppression factor S 2 /M 2 as compared with 
the twist 2 term, the overall numerical coefficients in Eqs. (|36|) and (|37j ) are of the same 
order at M 2 ~ 1 GeV 2 . 

The LCSR approach allows one to estimate the theoretical uncertainty of the predicted 
form factor. We did it in the following way. First of all, M 2 and S 2 were varied within 
allowed intervals. Furthermore, in order to investigate the sensitivity to the choice of 
the renormalization scale our calculation was repeated at two fixed scales Q 2 and M 2 
adopting the variation of the results as a theoretical uncertainty. Finally, accounting for 
the missing twist >6 terms we assume that the absence of the latter introduces an additional 
uncertainty equal to ±F^(Q 2 ). All abovementioned variations of the LCSR prediction for 
F n (Q 2 ) are then added linearly which is a rather conservative approach. 

In Fig. |^ we plot F n (Q 2 ), calculated with the asymptotic pion DA and at M 2 = 1 GeV 2 . 
The resulting uncertainty of the form factor indicated in this figure is about ±(20 -j- 30)% 
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at Q 2 > lGeV 2 . At Q 2 < 1 GeV 2 the uncertainty grows revealing the inapplicability of 
LCSR for small momentum transfers. In the region 2.0 < Q 2 < 10 GeV 2 our prediction 
for the pion e.m. form factor with the asymptotic pion DA can be fitted to the following 
simple formula: 

QW> = (0.0735 + 0.2016) + ^±^!2 _ 0.8496* 1.20K (ag) 

(all numbers in GeV 2 ), where the correlated intervals take into account the total uncer- 
tainty. 

The accuracy of the LCSR prediction (|32|) can be improved further by including various 
higher-twist corrections (due to twist-4 multiparticle and twist 6 DA) which were not yet 
analyzed. However, the smallness of the factorizable twist 6 term indicates that these 
effects are most probably numerically unimportant. In addition, it is desirable to improve 
also the perturbative expansion of the correlation function calculating the 0(a s ) term of 
twist 4 and the 0(a 2 ) term of twist 2. An attempt to account for the latter was made in 
Ref. [0 by matching LCSR to the NLO perturbative calculation |20 . 



Finally, in Fig. y we compare our numerical prediction with the recent accurate data on 
F-niQ 2 ) obtained from the pion electroproduction at Jefferson Lab [Tj| at Q 2 = 0.6 4- 1.65 



GeV 2 , in the region which only partly overlaps with the LCSR validity region Q 2 > 1 
GeV 2 . We find that within theoretical uncertainties and experimental errors the form 
factor calculated with the asymptotic pion DA ip n (u) is consistent with data |. 

With Fn^Q 2 ) accurately measured at the whole region Q 2 = 1 4 10 GeV 2 it should in 
principle be possible to constrain/fit the nonasymptotic part of <fn{u) determined by the 
coefficients a n in the expansion over Gegenbauer polynomials (see App. |BJ). Taking into 
account the complete expansion one obtains 

oo 

F^Q 2 ) = [F^{Q 2 )} as + J2 «2n(/io)/ 2 n(Q 2 ,/i,/io) , (39) 

n=l 

where the first term on the r.h.s. is the form factor calculated with the asymptotic DA 
and in the sum each am is multiplied by a calculable function: 

(°) la 

j2n{Q , V, Ho) = 6 I - j—^ I du uuCm [u - u) exp I ) + " ' ^ ^ 

In the above, /xq ~ 1 GeV is a certain low scale, and the anomalous dimensions 7„ of 
the renormalization factors are given in Appendix [BJ. For brevity, the 0(« s )-correction to 
Eq. fllOD is denoted by ellipses. 

A direct fit of all a n from Eq. (^) is of course not a realistic task. In fact, using the 
arguments of the conformal partial wave expansion, one expects that the coefficients are 
decreasing with n, a 2n+ 2 < &2n- Based on these arguments, the form of <p n (u) usually 

5 LCSR predictions also agree with older measurements of F 7r {Q 2 ) at Q 2 = 1 -f- 6 GeV 2 which however 
have large experimental errors. 
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Figure 4: Graphical illustration to Eq. ([3{\). The pion e.m. form factor obtained from 
LCSR with the asymptotic pion DA (solid line) and the coefficients at a2(lGeV) (long- 
dashed) and a^lGeV) (short- dashed). 



discussed in the literature involves one (02) or two (02, 04) nonzero coefficients neglecting 
the rest. Having adopted a certain simple ansatz for ip n (u) one is then able to constrain 
or even fit the coefficients from Eq. (|3"9"D. However, the current data |T8[ are sufficient 



to constrain only the simplest ansatz with a single nonzero coefficient 02- This can be 
seen from Fig. [| where /2 and f± in Eq. Q4TD are plotted in comparison with F£ s . Due to 
different signs of /2 and f± at Q 2 < 3 GeV 2 it is difficult to distinguish the form of ip n 
with <22 7^ from the one where both 02, a\ 7^ 0. E.g. one obtains equally good fits to the 
experimental data shown in Fig. (|J) with a2(lGeV) = 0.05, a^lGeV) = —0.30 as with 
ci2(lGeV) = 0.25, a 4 (lGeV) = 0. If we impose that all a n>2 = in Eq. (j3~9l) the coefficient 
<22 can be fitted to the following interval consistent with zero: 

a 2 (lGeV) = 0.24 ±0.14 ±0.08, (41) 

where the first error reflects our estimated theoretical uncertainty whereas the second one 
corresponds to the experimental errors. One needs data at larger Q 2 to resolve more 
complicated patterns of nonasymptotic coefficients in (p^. 
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5 The kaon electromagnetic form factor 



The LCSR for the charged kaon e.m. form factor can be easily obtained from the correlation 
function ([]), substituting P = K + , e\ = e u = +2/3, e2 = e s = —1/3, m qi = 0, m q2 = m s . 
We will systematically retain all 0(m s ) ~ 0(m 2 K ) effects which are numerous, in general. 
At the purely kinematical level one has to account for p 2 = m 2 K in the correlation function. 
Furthermore, the s-quark propagator produces a chirally non-invariant part proportional 
to m s which brings the twist 3 contribution into the game (the m 2 in the denominator 
of the quark propagator is neglected, being a higher-order effect). Finally, in the light- 
cone DA there are SU(3)-flavour symmetry (SU(3)fi) violating corrections of three types. 
First, the normalization factors, determining the quark- ant iquark vacuum-to kaon matrix 
elements in the local limit x = differ from the corresponding factors for the pionic 
matrix elements, e.g. fx ^ fn- Secondly, the nonasymptotic parts of the kaon DA are 
asymmetric with respect to the interchange of quark and antiquark fields, with a larger 
average momentum fraction of the strange quark. At the twist 2 level this effect manifests 
itself in the nonvanishing odd coefficients of the Gegenbauer expansion (|B7p: af , af , ... ~ 
m s 7^ 0. For the higher twist DA the SU(3)fi violating asymmetries were not studied yet, 
and in our numerical calculation we will neglect them. On the other hand, we will take into 
account the so-called meson-mass corrections to the twist 4 DA investigated and worked 
out in Ref. [JCT]]. These effects include the mixing of nonasymptotic parts of twist 2,3 and 
4 DA beyond the chiral limit. The corresponding expressions are presented in App. [B[ 

Apart from SU(3)fi violating corrections listed above the derivation of the LCSR for 
the kaon e.m. form factor repeats the procedure for the pion outlined in the previous 
section. The result reads: 

Fk(Q 2 ) = 4 2) (Q 2 ) + F^\Q 2 ) + 4 3) (Q 2 ) + 4 4) (Q 2 ) + 4 6) (Q 2 ) , (42) 
where the twist 2 contribution is 

2 1 A / uQ 2 um 



4 2) (Q 2 ) = j du (j<p K {w) + exp (-^ + ^fj ■ (43) 

The DA ip K (u,fi) is defined as in Eq. (|B1|), with q x = u and q 2 = s, so that u is the 
momentum fraction of the s quark. In the above, the lower limit u$ is related to the 
duality threshold in the kaon channel Sq by the equation = Uq(Q 2 /uq + m 2 K ) which 
should be solved with 0(m 2 K ) accuracy. The 0(a s ) correction to the twist 2 contribution 
has been calculated in Ref. in the chiral limit. To obtain F^' aa ^ we replace (p n by ifx 
in the expression for F^ ,as > given in App. 0. In addition, there are SU(3)fi violating 
corrections in the hard amplitude. Note that the first-order in m s corrections are absent 
due to chirality. Nevertheless, indirectly, 0(m s ) contributions will appear due to purely 
kinematical terms 0(p 2 = m 2 K ). To obtain these terms one has to recalculate the 0(a s ) 
diagrams retaining p 2 ^ 0, which is beyond the task of this paper. 
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The twist 3 and 4 terms in Eq. ( |42| ) can be cast in the same form as the twist 4 
contribution to the pion form factor: 



uQ 2 um 2 



'o 



uM 2 



+ 



A' 



{Q 2 -s«)m 2 K \ J3 , 



uM 2 M 2 



^Q2 + S K (Q2 + S K )3 



where 
and 



( 3 ) 2m s f 3K j , n , , 



(44) 



(45) 



- - ( 4 [2g 1K {u) + ^ia-(m)] - 4 [2G 2A -(w) + GW(«)] - 2u [2g 2K {u) - g 2K (u)} 



+ / Von 



1 - 2v) (2[2Lp ±K (oti) - (p± K {oti)] + 2ip\\ K (ai) - ^>\\k{ol^) 
-2 ^Kpucipti) + ^jjf(ai)] _ 2 ^||A-(ai) + <P\\k(oh 



.(46) 



terms are taken into account in accordance 



In the surface term in Eq. ( f44|) only 0(m 
with our approximation. Since the twist 3 contribution to the correlation function is 
proportional to m s it is consistent to use the SU(3)fi limit of the twist 3 DA, in particular 
¥>3K = (f3ir in Eq. (|4"5|). Finally, for simplicity we adopt SU(3)fi limit for the numerically 
small twist 6 factorizable term. 

In Fig. |5| we plot the kaon e.m. form factor calculated with the following choice of 
parameters: 

a) the twist 2 DA is taken with af(l GeV) = —0.17 (as estimated from the 2-point sum 
rule in Ref. |19[) and neglecting all higher Gegenbauer coefficients (in particular a 2 = 0), 
that is, maximally close to the asymptotic regime; 



b) 



1.2 GeV 2 is determined from the two-point QCD sum rules for fx [21]. Im 



portantly, Sq is larger than Sq reflecting heavier states in the kaon channel; 

c) for the strange quark mass we adopt an interval m s (lGeV) = 150 ± 50 MeV; 

d) the parameters of twist 3,4 kaon DA taken in the SU(3)fi limitQ: /^(lGeV) = 
/W(lGeV) = 0.0035 GeV 2 , 5 2 K (lGeV) = 5 2 (lGeV) = 0.17 ±0.05 GeV 2 (see Eq. (g§), 
u 3K {lGeV) = c^(lGeV) = -2.88, e 4 tf(lGeV) = e^(lGeV) = 0.5 0, [T|. The typical 



accuracy of all parameters except 5 2 is about 50%. We also use fx = 1.22/^. 

6 As noted above this is only consistent for the twist 3 part of the sum rule. The accuracy of the 
twist 4 part can be improved further if one determines the nonperturbative parameters entering the kaon 
twist 3,4 DA from the corresponding two-point sum rules in the kaon channel taking into account SU (3) /; 
violation, a task for future work. 
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Figure 5: The LCSR prediction for the charged kaon e.m. form factor Fx (Q 2 ) (solid line) 
in comparison with the pion form factor obtained with the asymptotic pion DA (dashed 
line) at M 2 = 1 GeV 2 . 



Comparing the LCSR prediction for the pion and kaon form factors calculated at M 2 = 
1 GeV 2 we observe a noticeable SU(3)fi violating difference. The ratio Fx(Q 2 )/F n (Q 2 ) 
approaches 1.5 at Q 2 ~ 10 GeV 2 . A closer look at Eq. ([43"D reveals that this difference 
originates from an interplay of two opposite effects. The SU(3) ji asymmetry due to a\ ^ 
in (Pk{u) tends to suppress the kaon form factor because in the larger contribution 
(corresponding to the u-quark interacting with the virtual photon) <Pk(u) < ¥n(u) in the 
end-point integration region Uq < u < 1. On the other hand, the fact that the duality 
threshold for the kaon is higher, sqk > So-k, implies that the end-point region for the kaon 
form factor is itself larger, thereby increasing Fk- Numerically, the latter effect turns out 
to be more important. Interestingly, the twist 3 contribution which is entirely an SU(3)fi 
violating effect is negligibly small, so that the m s uncertainty is unimportant. Note that in 
the ratio of kaon and pion form factors some theoretical uncertainties (e.g., due to M 2 - and 
scale-dependence) cancel leaving the major uncertainty in the Gegenbauer coefficient a\. 
The unaccounted SU(3)fi violating effects in higher twists can presumably be neglected 
within the present accuracy. If the kaon e.m. form factor is measured one would then be 
able to constrain/fit a\. 

Our final comment in this section concerns the neutral kaon e.m. form factor. It 
can be easily calculated from the same LCSR (^2j) if one replaces w-quark by (i-quark in 
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the initial correlation function, having in mind that DA of K° and K + are equal due to 
isospin symmetry. In particular, the leading twist 2 contribution is obtained replacing the 
u-quark charge 2/3 in Eq. (ff3|) by the d-quark charge -1/3. As a result F K o{Q 2 ) is a pure 
SU(3)fi violating effect proportional to the integral over the difference (px(u) — ipx(u). 
The numerical result is small : Q 2 F K o(Q 2 ) = 0.05 - 0.09 GeV 2 , (at 1 < Q 2 < 10 GeV 2 ) 
implying that the measurement of this form factor is a difficult task. 

We conclude that the LCSR method allows to systematically account for SU (3) break- 
ing effects in the kaon form factors, and that these effects revealed by the ratio of K + and 
7r + form factors are predicted to be quite noticeable. 

6 The K —> 7T form factor 

As a final application of LCSR in this paper we consider the K — > 7r form factor 
defined as 

<tt-(P " q) I si,u | K°(p)} = 2f^(q 2 )p, - (J+tf) - f^(q 2 ))q, . (47) 

As explained in Sect. 2 this form factor can be calculated from the correlation function (|I|) 
in two different ways, either from the vacuum-to-pion or from the vacuum-to-kaon corre- 
lation functions. Both calculations are valid only at sufficiently large spacelike momentum 
transfer Q 2 > 1 GeV 2 , whereas the form factor f£ n is measurable only at timelike momenta 
where the LCSR method is not applicable, e.g., at < q 2 < (mjf — m n ) 2 in the K e3 decays, 
or at {rriK + rn n ) 2 < q 2 < m 2 in r — > Kttu decays. 

Nevertheless, one is able to use the fact that LCSR obtained for two different settings 
yield one and the same physical parameter and derive useful constraints on the light-cone 
DA of the pion and kaon involved in both sum rules. To explain the idea we explicitly 
write down the LCSR obtained from the vacuum-to-pion correlator: 

&w2) ^/>"^" )exp (-S + S) + -- < 48 > 

Using instead the vacuum-to-kaon correlator, one gets 

fL{Q2) ^ T, C, du ^ "> exp ("S " 9? ) + - • < 49 > 

In the region 1 < Q 2 < 3 GeV 2 and at M 2 ~ 1 GeV 2 both sum rules are valid and the 
higher twist and 0(a s ) contributions denoted by ellipses are small, so that we may neglect 
them for the sake of simplicity. Equating fl48|) and (|49|) in this region one may constrain the 
nonasymptotic coefficients. In particular, the rate of SU(3)fi breaking asymmetry in the 
kaon DA ifx can be estimated if the pion DA <p n is determined with a sufficient accuracy. 

As a numerical illustration, in Fig|] we compare the r.h.s. of Eqs. and ( |I9"D 
calculated respectively with the asymptotic pion DA and with the kaon DA adopted in our 
calculation of Fk in the previous section, that is with ai(lGeV) = —0.17 and all a n> \ = 0. 
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g 2 [GeV z 



Figure 6: The K — > tx form factor at spacelike region calculated with twist 2 accuracy: 
from Eq. with the asymptotic pion DA (solid) and from Eq. ( $Q ) with the kaon DA 
including the SU(3)fi violating correction ~ a\ (dashed) and at a\ — (dash- dotted). 

The resulting agreement of two different sum rules is nontrivial and ensures confidence in 
the whole procedure and especially in the choice of duality thresholds in both pion and 
kaon channels. Note that the agreement is violated if we put a x = 0. 



7 Conclusions 

In this paper, we have studied the correlation functions of light-quark currents used to 
derive the LCSR for the pion and kaon form factors. We have demonstrated that the Ward 
identities for these correlators yield relations between DA of twist 4, a new alternative to 
using the QCD equation of motions. On the phenomenological side, we have corrected the 
expression for the twist 4 contribution to the LCSR for the pion form factor. The form 
factor calculated with the purely asymptotic pion DA is generally consistent with the recent 
Jefferson Lab data. On the other hand, constraining the nonasymptotic part of the pion 
twist 2 DA in terms of separate Gegenbauer coefficients demands more data at intermediate 
momentum transfers, 1 < Q 2 < 10 GeV 2 and largely depends on the particular ansatz 
adopted for this DA. A recent study of a similar problem for the 7*7* — > vr° form factor 
can be found in Ref. . 



We have presented the first LCSR prediction for the kaon e.m. form factor and demon- 
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strated that within the sum rule approach the SU(3)fi violating difference between kaon 
and pion form factors is systematically calculable in powers of the strange quark mass. It 
has been shown that a useful complementary information concerning the kaon DA can be 
obtained from the comparison of two independent sum rules for the K — > 7r form factor. 
In general, our results support the point of view that LCSR for the pion and kaon form 
factors combined with sufficiently precise data on these form factors represent a very useful 
tool for probing the pion and kaon light-cone distribution amplitudes. 
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A Light-cone expansion of the quark propagator 

The expansion of the quark propagator with a nonzero mass m near the light-cone (x\ 

x 2 ) 



reads (see e.g. Ref. 



S(xi 1 x 2l m) 



-i{0\T{q{x 1 )q{x 2 )}\0) 



d 4 k 



-ik(xi—X2) 



ft + m 

yy i 2 



rrr 



x 



(2vr) 4 

1 ft + m 

2 {k 2 — m 2 ) 2 ^ v k 2 — m 2 



- Jdv G (iv (vx 1 + (l-v)x 2 ) 



(Al) 



with = g s Gl u {\ a /2), Tr(A a A b ) = 25 ab . At m = 0, after the integration over k this 
expression reduces to the propagator derived in Ref. [p4] , (see also Ref. 

i 



S{x l ,x 2 ,0) 



27r 2 [(xi — x 2 ) 2 ] 2 167r 2 (xi — x 2 ) 



dvG^ivxx + (l-v)x 2 ) 



X 



(A2) 



B Light-cone distribution amplitudes 



The light-cone DA of the pseudoscalar meson P = n, K are defined according to Refs. 10 



TTJ. The matrix element of the axial- vector bilocal operator is expanded around the light- 
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cone (x\ = x\ = (x\ — x 2 ) 2 = 0): 



1 



(0|g 2 (x 2 )7^5gi(^i)|Pb)) = fp due- iu ^- iu ^{i P , ( VP {u) + (x, - x 2 ) 2 g lP (u)) 







+ (fa - x 2 ), - V ^ x \_ X ^ ) 92p(u)} , (Bl) 

retaining the leading twist 2 DA (fp(u) and the twist 4 DA g\p{u) and g 2 p(u), where u is 
the light-cone momentum fraction of the quark q^. From the local limit of Eq. ( |B1| ) one 
has the following normalization conditions: 

/ du(fp(u) = l, / dug 2 p(u) = 0. (B2) 
Jo J 

The twist 3 quark-antiquark DA ip pP and <p a p and the quark-antiquark-gluon DA ip 3P are 
defined as follows: 



1 



(0|fe(z 2 )i T5 <&(a:i)|P(p)) = fpfipjdue- iupx ^ px ^ p p(u) , 





2 



1 92(^2)^7551(^1) I Pip)) = — — ( 1 , 2 



ifpfip L mp 
6 V Vp 



1 

< Jdue- iupxi - iapX2 (p aP (u) , (B3) 




where \ip = m P /{m qi + m q2 ), and 

(0\q 2 (x 2 )cr^'j 5 G a(3 (vx 1 + vx 2 )qi(x 1 ) \P(p)) = if3p[(PaPn9pv - PpPp9av) 

-(PaPu9 P(t ~ Pf3P»9«»)} J Va,. <p 3P (a h ^e-i^-ia^-ia^+vp^) ? (fi4) 

all these DA being normalized to unity. Furthermore, the quark-antiquark-gluon twist 4 
DA are defined by the following matrix elements: 

(0|g 2 (x 2 )7 M 7BG^(^i + vx 2 )qi{xi)\P{p)) = fp J ^^-^—^—3(^1+^2) 
X \P» _ <P\\p( a i) + K9„aPp ~ 9^PaW±p{at)\ , (B5) 

(OMx 2 h l ,iG a p(vx 1 +vx 2 )q 1 (x 1 )\P( P )) = f P Jva-e-^-^-iasivp^+vp^) 

r p a (xi ~X 2 )p ~pp{xi -£ 2 )a~ ( x . ( X ± \~ f \\ ma\ 

x\Pn p ^ _ x ^ f\\Py a i) + \9^JPp - 9 tl /3Pa)<P±p{a i )\ , (B6) 
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where G a p = \e a f3 P \G px and the following abbreviations are used: 

(jCi - X 2 ) a Pf3 + (Xl - X 2 )/3P a 



T>ai = da\da2da^5 {l — a\ — a2 — a^) and g a/ 3 = g a /3 



p(x! - X 2 ) 



The distribution amplitudes are constructed |10[ using the formalism of the conformal 
expansion. The most familiar example is the twist 2 DA |i j 



(fp(u, jj) = Quu 



n=l 



(B7) 



where the expansion goes in Gegenbauer polynomials C^ 2 , the first four polynomials being 



Cl /2 (x) 

c 3 3 /2 (x) 



3x , G> 



3/2, 



'1 - 5x 2 



-x(3 — 7x) , C. 



3/2 



X 



1 Fi 

^-(1-Ux 2 + 21x 4 



(B8) 



The scale-dependence is given in the leading order by 



(B9) 



where L(fi 2 , Hi) = ot s (H2) / ot s (Hi) i (3q = 11 — \Np and the anomalous dimensions are 



7<°> = C f 



? jfi: 

(n+l)(n + 2) 



VJfc=l 



(BIO) 



For the pion, the coefficients vanish at odd n in the isospin symmetry limit. 

The twist 3 and 4 DA have been derived in Ref. |TU| using QCD equations of motion 



and conformal expansion. In Ref. |ITJ the meson mass correction have been worked out. 
We present here the explicit expression for these DA to next-to-leading accuracy in the con- 



formal spin, including the meson mass corrections (as explained in more detail in |10|, [TT 



and using the original notations of Ref. [Kj. Note that SU(3)fi violating nonasymptotic 
corrections to these DA (analogous to an 7^ for are still missing and have to be 
worked out in future. 

The twist 3 DA of the pseudoscalar meson, to next-to-leading order in conformal spin 
read: 



ip p p(u) 



1+ 30 



+ 



3P 



5 mi 



fjipfp 2 Hp J 
J3PUJ3P 27 m 2 



PaPW 

y?3p(ai) 



fipfp 20 Hp 
6uul 1+^5 



l + 6<))C 4 1/2 (2u-l), 



(Bll) 



Hp / 1 1 x 

HpJp 10 



^(l + far))C 2 3 / 2 (2n-l)},(B12) 



360aiQ;2«3 ( 1 + ^^-{7a>3 — 3) 



(B13) 
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The nonperturbative parameter f^p is given by the matrix element which corresponds to 
the local limit in Eq. ( |B4|) . The second parameter u 3 p determining the nonasymptotic parts 
of twist 3 DA is defined with the following matrix element (up to higher twist corrections): 



3 3 
(0\q^^\7B[Dp,Gax]qi - -jd^a^xlhGaxqiWip)) = j^fap^spPaPpPn ■ 

The scale dependence of the twist 3 parameters is given by: 

_j_ J-(I£f +n ) 

Hp(h 2 ) = [L(fj, 2 , Mi)] 00 VpM > fapM = [L(^2,fii)V oy - 3 ' f 3 p(Vi) 

(/m**p)M = [l(/x 2 ,^)]^v 6 + 3 ; (Apwsp)^; 

Finally, the four twist 4 three-particle DA defined in Eq. (|B5|), (|B6| ) are: 

9 



(B14) 



(B15) 
(B16) 



120(5pe P - i^a 2 mp){ai - a 2 )aia 2 a 3 



30(«i — a 2 )a\ 



-I206paia 2 a 3 



o P 



+ 25 2 P ep(l - 2a 3 



9 

40 



a 2 mp{\ - a 3 ) 



- + e P (l - 3a 3 ) 



0\\p(°n) 

<Px(ati) = 30a 2 
normalized as 

J Vaiip\\p(ai) = J Vai(p± P (ai) = , -J T>oaip\\p(aii) = J Pa^ip 



5 2 \ 9 

-y + 25pep(l - 2a 3 ) J (l-a 3 ) +— a,fmp(a 2 + a 2 - 4«ia 2 ) 



,(B17) 



r) 2 

f- (B18) 



The corresponding two-particle DA have the following expressions: 

— 2mm) — co>3P (3 — mm(27 — 56mm)) 



i \ 5 2 2 -2 , j /3^ m P ronci 



mp 



320 



5(25 - 29mm) - Y2al{\ - 5mm(19 - 52mm) 



>MM 



+ i(5pe-^a^mp) [2m 3 (10-15m + 6m 2 ) hiM + 2M 3 (10- 15m+6m 2 ) lnM+MM(2 + 13MM)](B19) 



92PW 



^p + ml (l + ^r(l-7MM))-^J(10-^ 3 p(l-TMM)) 



The nonperturbative parameter Sp is defined as 

(0\<l2G a ^ a qi\P(p)) = -iSpfpP» 



mm(m-m) .(B20) 



(B21) 
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with the scale-dependence: 



4(/i 2 ) = [L( M2)Ml )] ^ 5 2 P (m) 



(B22) 



whereas the second parameter ep determining the nonasymptotic corrections has the fol- 
lowing definition in terms of a local matrix element (up to twist 5 corrections) O ] : 

(0|g 2 [-D^, <5^]7 € gi - ^<9 M g 2 G^7 5 gi|P(p)) = —^fp5pe P (p^p v - ^mpg^ . (B23) 
The corresponding scale dependence is : 



(<5f>e P )Gu 2 ) = [L0*2,^i)] 3/3 ° (5pe P )(vi). 



(B24) 



C Radiative Corrections to the Twist 2 Pion Form 
Factor 

Here, for completeness, we present the formula for the radiative correction to the twist 2 
part of the sum rule for the pion form factor obtained in Ref . [[7] : 



i 

F (2, *.)(Q2 )= J du ^ u ^) 



Q(U~ Wo^softO, M 2 , S ) + 6(W - W^hardK M 2 , S ) 



(CI) 



where 



J" soft (M, M 2 , s 
a 



4vrM eXP (^) 



O 2 

9 + — + 3 In ^- + 3 In — - In 2 - In 2 



~3 



Q 2 



«g 2 



+ 



uQ 2 /u 



dsQ 2 e~ s / M2 
u{Q 2 + sf 



-P 



5s + I 1 + 2 In — £ I + 2 I — + s ) In — 



Q 2 



uQ 2 

U\l 2 



it 



2Q 2 /Q 2 + s 2M 2 + Q 2 + s, -p\ -p 
H 1 ttt; m I In 



m y s 

uQ 2 /u 



M 2 



+ 



rfs Q 2 e~*/ M 
ww(Q 2 + s) 3 



2« Q -s + sln 



s y /r 
s \ 



Q 2 + 5s + 2(Q 2 - s) In 



M 2 



u 



In — — In 

/I 2 M 



(C2) 



and 



J"hard(«,M ,Sq) 



24 



ds Q 2 e 



2 -s/M 2 



u(Q 2 + a) 



2[Q 2 -s + sln—) + - -Q 2 + 5s + 2(Q 2 



u 



s) In ^ 



s(Q 2 + s) 
M 2 



-3 + 2 In 



In 



UqUq /M 2 

uu 



2 In 



So 



3 In 



Po 
V 2 



(G3) 



Here p = uQ 2 — us and po = U Q 2 ~ usq = (1 — u/uq)Q 2 . 
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